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ABSTRACT 



Context. The single degenerate model is the most widely accepted progenitor model of type la supernovae (SNe la), in which a carbon-oxygen 
white dwarf (CO WD) accretes hydrogen-rich material from a main sequence or a slightly evolved star (WD +MS) or from a red giant star (WD 
+ RG), to increase its mass and explodes when approaching the Chandrasekhar mass. The explosion ejecta may impact the envelope of and 
strip off some hydrogen-rich material from the companion. The stripped-off hydrogen-rich material may manifest itself by means of a hydrogen 
. line in the nebular spectra of SNe la. However, no hydrogen line is detected in the nebular spectra. 
Aims. We compute the remaining amounts of hydrogen in red giant donors to see whether the conflict between theory and observations can be 
^ overcome. 

$_( Methods. By considering the mass-stripping effect from an optically thick wind and the effect of thermally unstable disk, we systematically 
ly-j carried out binary evolution calculation for WD + MS and WD + RG systems. 

, Results. Here, we focus on the evolution of WD -I- RG systems. We found that some donor stars at the time of the supernova explosion contain 
' little hydrogen-rich material on top of the helium core (as low as 0.017 Mq), which is smaller than the upper limit to the amount derived from 
' observations of material stripped-off by explosion ejecta. Thus, no hydrogen line is expected in the nebular spectra of these SN la. We also 
^ derive the distributions of the envelope mass and the core mass of the companions from WD -I- RG channel at the moment of a supernova 
' explosion by adopting a binary population synthesis approach. We rarely find a RG companion with a very low-mass envelope. Furthermore, 
our models imply that the remnant of the WD + RG channel emerging after the supernova explosion is a single low-mass white dwarf (0.15 
\0 Me - 0.30 Mo). 

Conclusions. The absence of a hydrogen line in nebular spectra of SNe la provides support to the proposal that the WD + RG system is the 
progenitor of SNe la. 

o 
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:'"! ' 1. Introduction from its companion, and explodes when its mass approaches 

rN ' the Chandrasekhar mass limit. The companion may be a 

H • Although type la supernovae (SNe la) are clearly important in ^^^^.^^^^^^^^ ,t^, (WD+MS) or a red-giant star (WD+RG) 

. y . determinmg cosmological parameters, e.g., Qm and (Riess (yungelson et al. mE Li et al. 7997, Hachisu et al. [199911 
et al. mm Perlmutter et al. [M, the progenitor systems of ^^^^^^ 2003, Langer et al. 2000] Han & 

SNe la have not yet been confidently identified (HiUebrandt p^dsiadlowski 2004; Chen & Li[2007 2009; Han 20083 Meng, 

& Niemeyer 2000; Leibundgut 2000). It is widely believed ^^^^^ ^ j^^^ 2009; Meng & Yang 2010a; Lu|2009i Wang et 
that a SN la is produced by the thermonuclear runaway of a 2009a 2009b 2010). The SD model has also been verified 

carbon-oxygen white dwarf (CO WD) in a binary system. The , , ^- / \/t o v i ^nmut a 1. 

■"^ ^ ■' ■' by many observations (see Meng & Yang|2010b|l. An alterna- 

CO WD accretes material from its companion to increase its ^. ■ »u n 1 u ui j ui j . /rM-,\ 

^ tive IS the theoretically less probable double degenerate (DD) 

mass. When its mass reaches its maximum stable mass, it ex- j 1 /tu o . 1 mo a m7 uu- 1 ino/i> • u- u 

model (Iben & Tutukov 1984; Webbink 1984), in which a sys- 

plodes as a thermonuclear runaway and almost half of the WD , r^r^wrr^ 1 u» 1 1 » u 

^ -' tem of two CO WDs loses orbital angular momentum by means 

mass is converted into radioactive nickel-56 (Branch 2004). f 1 j- -.• j c n t-u 

^ of gravitational wave radiation and finally merges. The merger 
Two basic scenarios have been presented. One is the single , 1 j vc . . 1 c . j 

^ ° remnant may explode if the total mass of the system exceeds 

degenerate (SD) model, which is widely accepted (Whelan & r^, , , , i- w »u ■ u tt h u j» 

° I 1 ^ ' ' J f \ lYie Chandrasekhar mass limit (see the reviews by HiUebrandt 

Iben [1973J Nomoto et al. i984 L In this model, a CO WD in- ^ Niemeyer|2000!and LeibundgutMl). 

creases its mass by accreting hydrogen- or helium-rich matter 

In the single degenerate model, the companion persists af- 

Send ojfprint requests to: X. Meng ter the supernova explosion. The supernova ejecta collides with 
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the envelope of and strips some hydrogen-rich material from 
the surface of the companion (Marietta et al. l2000l Meng, Chen 
& Hanl200T; Pakmor et al. 2008 ). The stripped-off hydrogen- 
rich material may manifest itself by means of narrow H„ emis- 
sion or absorption lines in later-time spectra of SNe la (Chugai 
11986 ', Filippenko 1997i ). The amount of the stripped-off mate- 
rial determines whether or not the narrow hydrogen line can 
be observed. Marietta et al. ( I2000I I ran several high-resolution 
two-dimensional numerical simulations of the collision be- 
tween the ejecta and the companion. They claimed that about 
0.15 - 0.17 Mq of hydrogen-rich material is stripped from a 
MS or a subgiant (SG) companion and 0.5 Mq from red giant 
(RG) star Meng, Chen & Han (2007) used a simple analytic 
method but a more physical companion model than that used 
in Marietta et al. (f2000) to simulate the interaction between 
SNe la ejecta and companions, and found that the minimum 
mass of the stripped-off material from a MS or SG star is 0.035 
Mq. However, the results of Marietta et al. (I2000I I and Meng, 
Chen & Han ( 120071 ) did not include confirmation by observa- 
tions, i.e., no hydrogen line was detected in nebular spectra of 
some SNe la and the upper mass limit to the stripped-off mate- 
rial was set to be 0.02 Mq (Mattila et al. l2005l Leonard |2007 ). 
Pakmor et al. (120081 ) used a more physical companion model 
than and a similar numerical simulation to that of Marietta et 
al. ( 120001 ) to recalculate the interaction between the supernova 
ejecta and companion. They found results similar to those of 
Marietta et al. (120001 ) . In certain circumstances, they claimed 
that these results agree with observations, and hence that the- 
ory does not conflict with observations. However, the special 
conditions envisaged by Pakmor et al. (2008) appear to be un- 
realistic according to detailed binary population synthesis re- 
sults (Meng & Yang 120 lOal . Based on the results of Pakmor 
et al. (120081 ). the amount of stripped-off material is between 
0.06 Mq and 0.16 Mq, which is consistent with the discovery 
of Marietta et al. (2000) and Meng, Chen & Han (120071 ) (see 
also Meng & Yang ,2010ai ). The results of Pakmor et al. ( 12008 1) 
therefore do not resolve the confliction between theory and ob- 
servations. Justham et al. f2009 ) proposed that the rotational 
effect of a CO WD may prevent its thermonuclear runaway un- 
til the accretion phase has ended, which could produce a RG 
companion with a low-mass envelope and reconcile theory and 
observations. We also suggest that failure to detect a hydrogen 
line in nebular spectra of some SNe la may imply that the WD 
+ RG channel is a means of producing SNe la. The amount of 
hydrogen-rich material obtained by Marietta et al. ( 2000 ) by 
means of WD + RG channel is higher than observed (0.5 Mq), 
which may be due to the simplistic RG model used by Marietta 
et al. aOOOi . 

In Sect, m we describe our binary evolution model. We 
present our evolutionary and binary population synthesis re- 
sults in Sect. [3] and our discussions and conclusions in Sect. 

El 

2. Model and physics inputs 

Meng & Yang (1201 Obi ) developed a comprehensive progenitor 
model for SNe la. In the model, the mass-stripping effect by op- 
tically thick wind (Hachisu et al. 1 19961 ) and the effect of a ther- 
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Fig. 1. An example of binary evolution calculations. The evo- 
lutionary track of the donor star in HRD is shown as solid 
curve and the evolution of orbital period is shown as dot-dashed 
curve. Dotted vertical line and asterisks indicate the position 
where the WD is assumed to explode as a SN la. 



mally unstable disk were included (Hachisu, Kato & Nomoto 
120081 Xu & Li 120091 ). The prescription of Hachisu et al. 
( |1999al ) for WDs accreting hydrogen-rich material from their 
companions was applied to calculate the WD mass growth. 
The optically thick wind and the material stripped-off by the 
wind were assumed to remove the specific angular momen- 
tum of WD and its companion, respectively. In Meng & Yang 
(1201 Obi ), both the WD + MS and WD + RG scenarios are con- 
sidered, i.e., Roche lobe overflow (RLOF) begins at either the 
MS or RG stage. After the RLOF, WD accretes hydrogen-rich 
material from the donor and increases its mass smoothly. When 
the mass of the WD reaches 1.378 Mq, the WD is assumed to 
explode as a SN la. The CO WDs may explode at the opti- 
cally thick wind phase or after the optically thick wind while 
in either the stable or unstable hydrogen-burning phase or the 
dwarf nova phase. They considered more than 1600 different 
WD close binary evolution and obtained a dense model grid 
for SNe la. Based on the comprehensive model, Meng & Yang 
(1201 Obi ) derived a Galactic birth rate of SNe la that is compara- 
ble to that inferred from observations. In the WD + MS chan- 
nel obtained in Meng & Yang ( 1201 Obi ), the companion models 
are similar to those in Meng, Chen & Han (^009 ). The masses 
stripped-off should then be similar to those found in Meng, 
Chen & Han (2007]), i.e., higher than O.O35M0. We therefore 
only considered the case of a WD + RG channel, which may 
only contribute to a fraction of all observed SNe la (see Meng 
& Yang |2010bl and Wang, Li & Han 2010). All our methods 
for calculating the binary evolution and the physics inputs into 
the binary evolution calculation are similar to those in Meng & 
Yang ( 1201 Obi ) (see Meng & Yang i2010b, for a detailed descrip- 
tion of our method). 
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3. Results 

3.1. Binary evolution 

In Fig.[T] we show the evolutionary track of the donor star in the 
Hertzsprung-Russel diagram (HRD, solid line) and the evolu- 
tion of the orbital period (dot-dashed line). The initial parame- 
ters of the binary system are also shown in the figure. The donor 
star evolves from the zero age main sequence (ZAMS). RLOF 
does not begin until the star enters the RG stage. The WD ac- 
cretes hydrogen-rich material from the donor and increases its 
mass to = 1.378Mq, where a SN la is assumed to oc- 

cur, and Mf* = 0.3125Mo and log(pS'^/day) = 1.8808. At this 
point, the donor star consists of a helium core of 0.2955 Mq 
and a very thin hydrogen envelope of 0.017 M©; (we refer to 
Fig. |2] which shows the evolution of donor mass, core mass, 
and hydrogen envelope mass, where the definition of the core 
is identical to that in Han, Podsiadlowski & Eggleton (119941) 
and Meng, Chen & Han (2008)). The envelope is so thin that 
the donor star has even left the RG and evolved to become 
a WD. Supernova ejecta impacts and strips off hydrogen-rich 
material from the envelope. Almost all the envelope material 
is stripped off (Marietta et al. l2000l l. The amount of the mass 
stripped-off from the companion should then be lower than 
0.017 Mq, which is lower than the upper mass limit of 0.02 
Mq obtained from observations (Mattila et al. i2005i: Leonard 
I2007I I. Therefore, no hydrogen line should be observed in the 
nebular spectra of this SN la, and our model can be reconciled 
with observations. This result may imply that the lack of de- 
tection of a hydrogen line in nebular spectra of some SNe la 
could be evidence that the WD + RG channel represents a way 
of producing SNe la. 

The RG model used by Marietta et al. (120001) consists of a 
helium core of 0.42 Mq and a thick envelope of 0.56 Mq. The 
RG model may deviate from reality because the model was not 
obtained from a detailed binary evolution calculation and mass 
transfer between CO WD and the RG star was not considered. 
If mass transfer had been considered, the structure of the com- 
panion would differ significantly from that used in Marietta et 
al. (120001 see Sects. ITU and 1431). 

3.2. Binary population synthesis 

To obtain the distributions of the envelope mass and the core 
mass of companion stars from WD + RG channel at the 
moment of SN explosion, we performed a detailed Monte 
Carlo simulation using Hurley's rapid binary evolution code 
(Hurley et al. 2000. ,2002i . In the simulation, if a binary sys- 
tem evolves to a WD + RG stage, and the system is located 
in the (log P\ M^) plane for SNe la at the onset of RLOF, we 
assume that an SN la is then produced. The envelope mass and 
the core mass of the WD binary at the moment of SN explo- 
sion are obtained by interpolation in the three-dimensional grid 
(M^p, M^, log P') obtained in Meng & Yang (I2010bl l. 

In the simulation, we follow the evolution of 40 million 
sample binaries. The evolutionary channel was described in 
Meng & Yang ( 12010b ). As for Meng & Yang (2010b), we 
adopted the following input to the simulation: (1) a constant 




Fig. 2. The evolution of secondary mass (sold line), core mass 
(dot-dashed line) and the hydrogen envelope mass (dashed 
line). Dotted vertical line indicates the position where the WD 
is expected to explode as a SN la. The zero point of time is set 
at the onset of RLOF. 
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Fig. 3. Snapshot distribution of the envelope mass Me and the 
core mass Mc of companion stars from WD + RG channel at 
current epoch with acE - 1 0. 

star formation rate (SFR)over the past 15 Gyr; (2) the initial 
mass function (IMF) of Miller & Scalo (1979); (3) the mass- 
ratio distribution is constant; (4) the distribution of separations 
is constant in log a for wide binaries, where a is the orbital sep- 
aration; (5) a circular orbit is assumed for all binaries; (6)the 
common envelope (CE) ejection efficiency q-ce, which denotes 
the fraction of the released orbital energy used to eject the CE, 
is set to be either 1.0 or 3.0 (see Meng & Yang (2010b) for 
details of the parameter input). 

In Figs |3] and |4] we show current-epoch-snapshot distribu- 
tions of the envelope mass and the core mass of companions 
produced by the WD + RG channel at the moment of super- 
nova explosion with different acjH In these figures, we can 
discern a clear trend, i.e., that core mass decreases in general 



' In the figures, we do not show the cases of WD + MS because: 
(i) the definition of the core for a MS star is incorrect; (ii) the interac- 
tion between supernova ejecta and the companion is complex, i.e., the 
stripped-off mass from the companion strongly depends on both the 
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Fig. 4. Similar to Fig.[3]but for uce - 3.0. 

with envelope mass. This is a natural consequence of stellar 
evolution. In addition, varying ckce may have a significant ef- 
fect on the distributions. For acE - 1.0, the envelope mass of 
companions is mainly between 0.1 Mq and 0.6 Mq, while it 
is lower than 0.4 Mq for q-ce - 3.0. This is because that for 
a given primordial binary system, high qtce is indicative of a 
longer orbital period after CE ejection. When the RG star fills 
its Roche lobe, the star then consists of a low-mass envelope 
and a high-mass core. As a result, the lower mass limit of the 
core mass for acE = 3.0 is higher than that for acE - 1-0, 
while the upper mass limit of the envelope mass for acE = 3.0 
is lower than that for acE - 1-0. The core mass is between 
0.15 Mq and 0.3 Mq. After the interaction between supernova 
ejecta and RG companion, the RG companion loses almost its 
entire envelope (96% - 98%) leaving only the core of the star 
(Marietta et al. l2000l l. This may be a channel to forming single 
low-mass white dwarf (Justham et al. l2009l l. 

However, there are only a few systems with Mg < O.O2M0 
in our sample, which means that the binary sequence consid- 
ered in this paper may only be able to explain some individual 
observations (Mattila et al. .2005. Leonard l2007l l. 

4. Discussions and conclusions 

4.1. Age 

We have found that if a SN la originates from the WD + RG 
channel, the donor star may almost become a helium WD and 
only some hydrogen-rich material remains on top of the he- 
lium WD (as little as 0.017 Mq), which means that the upper 
limit mass to the stripped-off from the companion by super- 
nova ejecta is 0.017 Mq. No hydrogen line is then expected in 
the nebular spectra of some SNe la and our results may explain 
the conflict between theory and observation, i.e, theory predicts 
that the stripped-off material should be greater than 0.035 Mq 
(Marie tta et al. l2000l Meng, Chen & Han [20021 Meng & Yang 
I2010bl l. while observations indicate that the upper mass limit 

evolutionary stage of the MS star at the onset of RLOF and the initial 
parameters of a WD -1- MS system (Meng, Chen & Han 2007); and 
(iii) the study of WD + MS is beyond the scope of this paper. 



of the stripped-off material is 0.02 Mq (Mattila et al. 120051 
Leonard I2OO7I 1. No hydrogen line has been detected in nebu- 
lar spectra of some SNe la may indicate that the progenitors 
of the observed SNe la are from WD + RG systems. If the 
observed SNe la (SN 2001el, 2005am, and 2005cl) were pro- 
duced by the WD -1- RG channel, they should originate in an 
old population. Unfortunately, there is no constraint on the age 
of the three SNe la. We checked the types of their host galaxies 
and found that apart from the host galaxy of SN 2005cl (MCG- 
01-39-003), which is a SO galaxy (Wang |2009l Bufano et al. 
[2009), they are both spiral galaxies, i.e. the host galaxy of SN 
2001el (NGC 1448) is a Scd galaxy (Wang et al.'2003', Wang 
et al. 2006 ) and the host galaxy of SN 2001am (NGC 281 1) is 
a Sa galaxy (Bufano et al. |2009] l. The progenitor of SN 2005cl 
may therefore belong to an old population, but we are unable to 
infer any information about the population of the other SNe la. 
However, we note that all three SNe la are located at the edge 
of their host galaxy and that SN 2005cl is even located in a tail 
extending from MCG-0 1-39-003. Is this phenomenon evidence 
of an old population? It is possible because halo stars in general 
belong to an old population. All three SNe la are located at the 
edge of their host galaxy may be an observational select effect 
because we are more likely to observe a SN la at the outskirts 
of a host galaxy rather than its inner part. This selection effect 
might increase the probability that a SN la is produced by the 
WD H-RG channel with a low-mass envelope being observed. 

4.2. WD + RG system 

Relative to that of the WD -1- MS system, the Galactic birth 
rate of the WD + RG channel is low (see Meng & Yang |20I0bl 
and Wang, Li & Han (1201 Oi l). However, since the Galactic birth 
rate of SNe la predicted by the model in Meng & Yang |20I0bl is 
lower than that inferred from observations, the WD + RG chan- 
nel should be carefully investigated because the progenitors of 
some SNe la (e.g. SN 2006X and SN 2007on) are possible WD 
4- RG systems (Patat et al. 12007a ; Voss & Nelemans |2008] |. 
In addition, some recurrent nova (belonging to WD + RG) are 
suggested to be the candidates of SNe la progenitors (Hachisu 
et al. 1 999b, Hachisu & Kato 2006 ; Hachisu et al. lIOOTl l. In this 
paper, we even proposed that the prevalence of the WD + RG 
channel is why no hydrogen line was detected in nebular spec- 
tra of some SNe la, although the probability of its occurrence 
is low. 

4.3. Interaction between supernova ejecta and 
companion 

Marietta et al. (2000) performed several high-resolution two- 
dimensional numerical simulations of the collision between the 
supernova ejecta and companion, and found that a red-giant 
donor will lose almost its entire envelope (96% - 98%) due to 
the impact leaving only the core of the star (^ 0.42Mq). The 
RG star used in Marietta et al. (120001) consists of a helium core 
of 0.42Mo and an envelope of O.56M0, which are not compara- 
ble to those obtained by our simulations (see Figs[3]and|4|. In 
addition, the radius of their RG model is 180 Rq, which corre- 
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spends to an orbital period of ~ 900 days. Too long to compare 
with our simulation (see Fig. [Hi, the orbital period leads to a 
lower envelope binding energy than produced by the model de- 
veloped in this paper since the binding energy of the envelope is 
determined mainly by the radius of the RG star (Meng, Chen & 
Han l2008l l. The envelope of the RG model used by Marietta et 
al. (I2OOOI 1 is then more likely to be stripped off and the amount 
of material stripped-off by the RG companion in Marietta et 
al. (120001) might be overestimated. A more detailed numerical 
simulation of the interaction between supernova ejecta and an 
RG companion should therefore be performed by a more phys- 
ical companion model than that in Marietta et al. (I2000I I. 

4.4. Alternative explanation of the lack of hydrogen 

The absence of hydrogen lines in the nebular spectra of SNe 
la may have other explanations. The RG companion with a 
small hydrogen-rich envelope may be the results of either 
a fine-tuning effect as suggested in this paper, or a physical 
process. For example, Justham et al. ( I2009I I suggested that 
a rotational effect of WD could prevent the thermonuclear 
runaway occurring until the accretion phase has ended, which 
could also produce a RG companion with a low-mass envelope. 
Rotation may also increase the probability of a SN la from 
WD +RG channel with low-mass envelope being observed. 
An alternative explanation of the lack of hydrogen is that the 
amount of stripped-off material might has been dramatically 
overestimated as discussed above. 

Based on the discussions above, further attempts to observe 
hydrogen lines in nebular spectra of SNe la are encouraged to 
check our suggestion. The WD + RG system may also be an 
origin of single low-mass white dwarfs. 
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